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In this paper we prove that the inequality introduced by Collins, Gisin, Linden, Massar and 
Popescu [11] is tight, or in other words, it is a facet of the convex polytope generated by all 
local-realistic joint probabilities of d outcomes. This means that this inequality is optimal. We 
also show that, for correlation functions generalized to deal with three-outcome measurements, the 
satisfyability of this inequality is a necessary and sufficient condition for the existence of a local- 
realistic model accounting for them. 
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I. INTRODUCTION 

o ■ 

Quantum Mechanics predicts that there are events in Nature exhibiting correlations which are not explainable in 
terms of local-realistic theories [1,2]. Those correlations can arise when measuring two or more separated systems 
which are in an entangled state. This fact is known as quantum nonlocality and has been experimentally proven 
up to some loopholes [3] . Knowing which multipartite quantum states violate local-realism is a very important open 
t— I ' problem. On the other hand, the sharing of non-local-realistic correlations between some parties is a useful resource for 
tasks like distributed computation [4] and secret communication [5] . Then, studying whether a given set of correlations 
are achievable with classical randomness or, contrary, they need of quantum entanglement, is an important issue. 

It was shown in [6] that the sets of local-realistic joint probability distributions are convex polytopes. There is 
one of these polytopes for each setting. By a setting we mean a fixed number of: parties, observables per party and 
number of outcomes per observable. These polytopes can be characterized by a finite number of linear inequalities 
, that the joint probabilities of the correlated events must satisfy. These inequalities correspond to the facets of the 
polytopes, here we call them tight Bell inequalities. There are algorithms that find all the facets of a polytope, but 
the time they need for doing so grows very fast as the setting becomes less simple. Thus, obtaining these inequalities 
is in general a very hard problem, and it only has been completely solved in simple two-outcome settings [8,9]. In 
the case where correlation functions instead of joint probabilities are used, the problem is also completely solved 
for the setting consisting of an arbitrary number of parties each possessing two dichotomic observables, in a very 
mathematically-elegant way [10]. In this paper we give a way for generalizing correlation functions to deal with 
d-outcome measurements. Using this instead of joint probabilities makes the problem of finding all the facets of the 
corresponding polytope numerically more feasible. We have done an algorithm that finds all the inequalities that 
characterize the set of local-realistic generalized correlation functions for any d. It turns out that for d = 2,3, all the 
inequalities are equivalent to the CGLMP-inequality. 

In the literature Bell inequalities exist which are not facets of correlation polytopes (non-tight Bell inequalities), 
■ they are half-spaces that contain the polytope but do not lie in its frontier. In this sense we say that, tight Bell 
inequalities are the optimal detectors of non-local-realistic correlations. For most of the Bell inequalities there is no 
proof of tightness or non-tightness. In this paper it is shown that the CGLMP-inequality [11] is a tight Bell inequality. 
In other words, the CGLMP-inequality is a facet of the polytope corresponding to the setting of two parties with two 
d-outcome observables per party. In fact, it is a family of equivalent facets of this polytope, but not all; thus, we do 
not completely solve the problem of characterizing all local-realistic joint probabilities for our setting, we just give a 
necessary condition for them. 

This result is also important because often, the degree of violation of an inequality is used as a measure of how 
non-local a quantum state is. But this does not make sense unless a tight Bell inequality is used, since otherwise, the 
measure has some inherent bias. 

Moreover, as it was shown in [15], the CGLMP-inequality (for d — 3) is not maximally violated by the maximally 
entangled state when orthogonal measurements are done. This result becomes very surprising after knowing that this 
inequality is tight. 

This paper is organized as follows. In section II we show that the set of all local-realistic correlations is a convex 
polytope, and we characterize it in terms of its generators. Section III contains the proof of tightness for the CGLMP- 
inequality. In section IV we simplify the problem in a way that makes possible to find numerically all the inequalities 
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for d = 3. Some conclusions are commented in section V. In the appendix, the two lemmas used in section II are 
proven. 

II. CHARACTERIZATION OF THE POLYTOPE 

In this section we justify the use of convex polytopes when studying local-realistic correlations. We also explain the 
dual description of a convex polytope: in terms of its generators on one side, and in terms of its facets on the other side. 
The setting that concerns us involves two parties: Alice and Bob. Alice can carry out two possible measurements, A\ 
or A2, and Bob can carry out B\ or B 2 . Each measurement has d possible outcomes: Ai, A 2} Bi, B 2 = 0, . . . , d — 1. 
In the remaining of the paper the letters characterizing the measurements are also going to be used as the variables 
denoting their results. For each of the four experimental settings (A\B\, A\B 2 , A 2 B\, A2B2) there is a joint probability 
distribution of the outcomes: 

P(A a = k, B b = s) k,s = 0,l,...d-l a, 6 =1,2. (1) 

We can arrange all those 4d 2 numbers in a column vector, P <E lZ 4d2 , in order to have a geometrical formulation of 
the problem. For each experimental setting its corresponding joint probability (1) is normalized, hence, P must fulfil 
the four constrains 

d-i 

P(A a = k,B b = s) = l, a,6=l,2. (2) 

k,s=0 

Here, we are only interested in correlations that can not be used for instantaneous communication between Alice and 
Bob. That is, marginal probabilities for one party are independent of the measurement chosen by the other party, 

d-l d-l 

Y J P{A a = k,B 1 =s)=Y J P{A a = k,B 2 = s) . (3) 

This constraint must hold for each observable, A\, A 2 , B\, B2, and for each outcome k = 0, . . . d— 1, but not all of them 
are linearly independent. In Lemma 1 (in the appendix) it is proven that the normalization (2) plus non-signaling (3) 
conditions form a system of 4d linearly independent equations, which will be used later. Then, all possible P vectors 
— that satisfy these constrains — belong to an affine space [12] of dimension 4d(d — 1). 

In what follows, we characterize all joint probability distributions P(A a =k,B\ = s) obtainable with local-realistic 
models. Let A label all possible outcomes of the measurements, that is, the d 4 possible values of (Ai,A 2 , Bi,B 2 ). As 
it is well known, all these models can be written as 

P(A a = k, B b = s) = J2 PA P(A a = k, A) P(B b = s, A) with p x > and = 1, (4) 

A A 

where P(A a = k, A) and P(B b = s,A) take only the values 0,1. In other words, P(A a — k,B b = s) is a convex 
combination of extreme product probabilities, P(A a = k,B b = s) = 8A a ,k §B b ,s, which when written as column 
vectors are denoted by Ga- In our geometric picture, 

P = Pa G A , (5) 

A 

is equivalent to say that P belongs to the convex hull (Cg) expanded by the set of vectors Ga- We say that each one 
of the d 4 vectors Ga is a generator of Cg- Because the number of generators is finite, Cg is a convex polytope. It 
is easy to see that the generators satisfy the normalization (2) and the non-signaling (3) conditions, and this implies 
that all the points belonging to the polytope also satisfy them. In Lemma 2 (in the appendix) it is proven that the 
affine hull [13] expanded by the polytope has dimension h = Ad{d — 1), just the maximum allowed by these conditions 
(2-3) as proved in Lemma 1. It is important to know the dimension of this space, h, for what follows. 
It is known that every convex polytope is characterized by a unique finite set of inequalities [14], called facets: 

PGC G ^ X r P< Xj j = 1,2... (6) 

Each of these facets fulfil the next two conditions: 
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Condition 1: Every one of the generators Ga must belong either to the half-space Xj • Ga < Xj or to the hyperplane 



Xj • G A = Xj. 

Condition 2 (tightness): Among the generators that belong to the hyperplane Xj - Ga = Xj there must be h which 
are affinely independent. Notice that an hyperplane of dimension h — 1 is completely characterized by h affincly 
independent points. In our particular case, the null vector does not belong to the polytope, hence the condition 
of existence of h affinely independent vectors is equivalent to the existence of h linearly independent vectors. 

All these inequalities can be grouped in families of equivalent inequalities. Two inequalities are equivalent if we can 
transform one into the other by composing the following symmetry transformations: 

Party exchange: P(A a = k,B b = s) — ► P(A b = s,B a = fc) or (A a B b ) — > (A b B a ). 



Observable exchange: P(A a = k,B b = s) — ► 

Relabeling of outcomes: P(A a = k,B b = s) 
permutation of the set {0, 1 ... d — 1}. 



P{A a = k,B b = s) or (A a B b ) — > (A a B b ), where 1 = 2 and 2 = 1. 
— ► P(A a = <r k ,B b = s) or (A a B b ) — ► -(A a B b ), where a k is a 



Therefore, the obtention of an inequality automatically yields all its family of equivalent ones by transforming it. 

Some of these inequalities follow from the fact that the joint probabilities satisfy the normalization and non-signaling 
conditions, thus quantum joint probabilities never violate them. We call them trivial inequalities. The non-trivial 
ones are the tight Bell inequalities and could be violated by quantum correlations. 

Non-tight Bell inequalities are the ones that fulfil Condition 1 but not Condition 2. They are worse detectors 
of non-local-realistic correlations because they could be modified in order to detect the same points as before plus 
additional ones. In the next section we show that the CGLMP-inequality is tight. 



III. THE CGLMP-INEQUALITY 

In this section we will show that the CGLMP-inequality satisfy condition 1 and condition 2 stated in the previous 
section. If we use the symbol = to denote equality modulus d, the CGLMP-inequality for any value of d can be 
written as 



MP) 



Ld/2-lJ 

£ i 1 

k=0 



2k 
d-l 

+P(A 2 



P(A! — Bi = k) - P{Ax - Si = -k - 1) + P{A 1 -B 2 = -k) - P(Ai -B 2 = k + l) 



Bt 



-k - 1) - P(A 2 -Bi = k) + P(A 2 -B 2 = k)- P{A 2 — B 2 = —k — 1) 



<2, (7) 



where [^J denotes the largest integer less or equal than x. Notice that for d — 2 (7) is equivalent to the famous 
CHSH-inequality [2]. The main result of this paper is that, for all d, the CGLMP-inequality is tight. In order to prove 
it, we need to show that both condition 1 and 2 are satisfied by (7). In [11] the authors showed that condition 1 is 
satisfied by the inequality. In what follows, we will reproduce their proof for the sake of completeness. 

Proof of Condition 1: Let us start by defining the four variables 

r = A 1 -B 1 +d, s= -At +B 2 + d , t= -A 2 + B 1 -l + d, u = A 2 - B 2 + d , (8) 

were d denotes a multiple of d (possibly different for each variable) that has to be added in order to have that 

d-l 



<r,s,t,u< 



(9) 



is satisfied. Recalling the definition of = and d, it is easy to see that this new variables fulfil the constrain 

r + s + t + u=-l . (10) 
Using these variables we have defined, the inequality (7) can be written as 



Ld/2-lj 



k=0 



2k 



P(r = k) - P(r = -k - 1) + P(s = -k) - P(s = k + 1) 
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+ P(t = -k - 1) - P(t = k) + P(u = k) - P(u = -k-l) 



(11) 



Where, because we have defined r, s, t, u such that (9) is satisfied, the = symbols in (7) can be substituted by simple 
equalities. After some algebra, the value of I d when applied to the generators, G\, yields 

J d (G A ) = /(r) + /(*) + /(*) + /(«) f{x) = { , (12) 

L d-i d-i ' x ^ u 

where the values of r,s,t,u are the ones corresponding to each A. Now, to prove Condition 1, we look at the four 
cases that appear when considering the different sings the variables r, s, t, u can have: 

Case 1. r, s, t, u are all positive. Then, (9) and (10) imply that r + s + t+ l = d — 1, which when inserted into (12) 
gives I d = 2. 

Case 2. Three of the variables r, s, t, u are positive and one is strictly negative. Then, (9) and (10) imply that either 
r + s + t + u = d — 1 which gives us Id = — 2/(d — 1), or, r + s + t + u = — 1 which give us Id = 2. 

Case 3. Two of the variables r,s,t,u are positive and the other two are strictly negative. Then, (9) and (10) imply 
that r + s + t + u = —1, giving Id = — 2/(d — 1). 

Case 4. One of the variables r,s,t,u is positive and the other three are strictly negative. Then, (9) and (10) imply 
that either r + s + t + u= —1 ovr + s + t + u= — d — 1 happens. Then, the two possible values for (12) are 
either I d = -2{d + l)/(d - 1) or I d = -2/{d - 1). 

Case 5. r, s, t, u are all strictly negative. Then, (9) and (10) imply r+s+t+u = —d—1 and then Id = — 2(d— 1). 

We have seen that for all the generators, G\, the value ^(Ga) is equal or less than 2, which completes the proof of 
condition 1. 

Until now, it has been proven that the CGLMP-inequality is a Bell inequality, that is, any set of joint probabilities 
that violates it is not inside the polytope. In what follows, we will see that this inequality (7) is a facet of the polytope, 
or in other words, it is as tight as it can be. 

Poof of Condition 2 (tightness): As we have seen in the proof of Condition 1, the generators contained in the 
hyperplane — that is, those that satisfy I d {G\) — 2 — are (i) the ones having the four variables r,s,t,u positive, 
and (ii) the ones having three variables positive and one strictly negative that satisfy r + s + t + u= —1. Now, we 
will demonstrate that among these generators there are h which are linearly independent. First notice that all the 
generators can be written as 

G(a u a 2 ,b u b 2 ) - |^i, Si) 8 L4i, B 2 ) ® \A 2 ,Bx) \A 2 ,B 2 ) , (13) 

where \A, B) stands for \A mod d) (g> \B mod d), and |n) is a c?-dimensional vector with a 1 in the n th component and 
0s in the rest 1 . Of course, every one of the four spaces in the direct sum has dimension d 2 . Using the variables defined 
in (8) we can also write (13) as 

\A 1 ,A 1 -r) 8 \A 1 ,A 1 + s) ®\Ax-r -t-l,A 1 -r) @ |A X +s + u,A 1 + s) . (14) 

Because we are only interested in the linear dependence properties of these vectors, we can apply to them a linear 
transformation that preserves orthogonality. Thus, (14) can be transformed into 

\Ax,r) 8 \A u s) |Ai - r,t) 8 |^i + s,u) , (15) 

by applying the permutation X)l _ fe= (IA k )(A A - k\ <g> \A, k)(A, A + k\ <g> \A, k)(A - k - 1,A\ ® \A,k)(A + k,A\). 
which is an orthogonal matrix. The construction of the set of Ad{d — 1) linearly independent vectors is done step 
by step. In each step, Ad vectors which are linearly independent among them and among the previously introduced 
vectors, are added. The following two examples show how to introduce a set of Ad vectors, for two particular cases, 
that will be used later for the general case. 



J To avoid confusion we remark that the order of the components in (13) is not the same than in (34), in Lemma 2. 
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Example 1: Imagine that only vectors of the form (15) having the values of r,s,t,u different from a have been 
introduced. Then, the following Ad vectors 



\A, a) © \A, 61) ®\A-a, b 2 ) ®\A + 61, 63) 
\A, 63) © \ A, a) © \A - 63, h) © I A + a, b 2 ) 
\A, b 2 ) © \A, 63) ®\A- 6 2 , a) © \A + 63 A) 
\A, 61) © |,4,62) © |A - 61,63) © \A + b 2 ,a) 



for A = 0,l...d-1, (16) 



where 61,62,63 + 1 a, can be added to the set, with the certainty that they are linearly independent form the 
previously introduced ones — since they did not have components like \A, a) in any of the four orthogonal 
subspaces — , and that they are also linearly independent among them, because \A, a) is, for each vector, in 
a different orthogonal subspace. Notice that argument is valid for all the values of 61,62,63 different from a, 
independently from the fact that they have appeared in a previously introduced vector, or not. 

Example 2: Imagine that only vectors of the form (15) having the values of r,s,t,u different from a have been 
introduced. Then, the following 4d vectors 



\A, a) © \A, a)®\A- a, 61) © \A + a, 6 2 ) 
\A, a) © \A, 61) © \A - a, a) © \A + 61, 6 2 ) 
\A, a) © \ A, 61) ®\A-a, b 2 ) ®\A + 61, a) 
\A, 61) © \A, a) © \A - 61, a) © \A + a, 6 2 ) 



for A = 0, 1 . . . d - 1 (17) 



can be added to the set, with the certainty that they are linearly independent form the previously introduced 
ones, for the same reason than in example 1. To see that they are also linearly independent among them, we 
can observe that for a fixed value of A, the projection of these four vectors onto the subspace \A, a)(A, a\ © 
I A, a)(A, a\ © \ A — a, a)(A — a, a\ © \ A + 61, a)(A + 61, a\ contains the minor 



110 
10 10 
10 1 

10 0, 



(18) 



which has a non-zero determinant. The argument is also valid for all the values of 61,62 different from a, 
independently from the fact that they have appeared in a previously introduced vector. 

Now, we proceed to construct the set of 4d(d—l) linearly independent vectors. As we said, we will do this operation 
in d — 1 steps. In each of these steps, 4cZ vectors are introduced in the same way than in examples 1 and 2, but, 
with different values for 0,61,62,63 in each step. To simplify this operation, we study separately four cases, namely: 
d = 4e, 4e + 1, 4e + 2 and 4e + 3, where e is an arbitrary integer. 

Case d = 4e. In this case the constrain (9) is — 2e < r, s, t, u < 2e — 1. For the first 2e — 1 steps, we introduce vectors 
of the form (15) with r, s,t,u> 0, for which (10) isr + s + f + M = 4e— 1. In each step, the scheme of example 
1 (16) is used with the following values for a, 61, 62, 63 of: 

a e-l e+1 e-2 e+2 e-3 e+3 e-k e+k 

61 e e — 1 e + 1 e — 2 e + 2 e — 3 e + k— 1 e — k 2e — 1 

6 2 ~ e ' e - 1 ' e'e-1' e'e-1'"' e 'e-l' - ' e 1 ' 

6 3 eeeeee e e e 

Notice that in each step — reading from left to right the top value of each column appears for the first time, 
hence the scheme of example 1 can be applied in every step. After this operation we have a set of 4e?(2e — 1) 
vectors. For the next 2e steps, we introduce vectors of the form (15) having three of the variables r,s,t,u 
positive and one strictly negative that fulfil r + s + t + u= —1. We also use the scheme of the first example but 
with the following values for a, 61, 62, 63 : 

-2 -3 -k -2e 

(20) 



Notice that as before, the top value of every column appears also for the first time. After this, we have 
constructed a set of M(d — 1) linearly independent vectors. This finishes the proof for the case d = 4e. 



a 


-1 


-2 


-3 


-k 


-2e 


61 





1 


2 


k - 1 


2e - 1 


62 


' 


' 


■• 


■• 





b 3 


















5 



Case d = 4e + 1. In this case (9) is — 2e < r, s, t, u < 2e. For the first 2e steps, we introduce vectors with r, s, t, u > 
and r + s + t + u = 4e in the following order: 







,y i 


e 4- 1 


e-1 










6 2 


e + 1 


e + 1 










a 


e-1 ' 


e 










a 


e - 1 


e 






a 


e + 2 


e - 2 


e + k 


e — fc 


2e 





h 


e-1 


e + 2 


e-k + 1 


e + k 


1 


2e 


b 2 


e - 1 ' 


e 


e-1 


e •■ 


■ e-1 ' 


e 


b 3 


e 


e 


e 


e 


e 


e 



(21) 



The first two steps have to be done following the second example, because each one has two times a value that 
appears for the first time. For all the rest, the top value appears for the first time, thus the first example can 
be applied as model. After this operation we have introduced 4d2e vectors in our set. For the next 2e steps, we 
introduce vectors having three of the variables r, s, t, u positive and one strictly negative and r + s + t + u = — 1. 
We also use the scheme of example 1 in the next steps: 



a 


-1 


-2 


-fc 


-2e 


bi 





1 


fc- 1 


2e- 1 


b 2 


' 


o •■ 


o •■ 





b 3 















(22) 



We have constructed a set of Ad(d — 1) linearly independent vectors, which finishes the proof for the case 
d = 4e + l. 

Case d = 4e + 2. In this case — (2e + 1) < r,s,t,u < 2e. For the first 2e steps, we introduce vectors which have 
f, s, t, u > and r + s + t + u = 4e+l, using the scheme of example 1, in the following order: 



a e + 1 


e - 1 


e + 2 


e - 2 


e + fc 


e — fc 


2e 





&i e 


e + 1 


e-1 


e + 2 


e-k + 1 


e + fc 


1 


2e 


b 2 e ' 


e + 1 ' 


e ' 


e + 1 


e 


e + 1 


■ e ' 


e + 1 


63 e 


e 


e 


e 


e 


e 


e 


e 



(23) 



For the next 2e + 1 steps, we introduce vectors having three of the variables r,s,t,u positive and one strictly 
negative and r + s + t + u = —1, using the scheme of example 1, in the following order: 



a 


-1 


-2 


-fc 


-(2e + l) 


bi 





1 


fc - 1 


2e 


62 


' 


•■ 


"' 





6.3 















(24) 



We have constructed a set of Ad(d — 1) linearly independent vectors, which finishes the proof for the case 
rf = 4e + 2. 

Case d = 4e + 3. In this case — (2e + 1) < r,s,t,u < 2e + 1. For the first 2e + 1 steps, we introduce vectors with 
r, s, t, u > and r + s + t + u = 4e + 2 in the following order: 









a 


e + 1 












bi 


e 












a 


e + 1 












b 2 


e 






a 


e - 1 


e + 2 


e - 2 


e + fc 


e — fc 


2e + 1 


h 


e + 1 


e-1 


e + 2 


e - fc + 1 


e + fc 





b 2 


e + 1 ' 


e + 1 ' 


e + 1 ' 


e + 1 ' 


e + 1 •' 


' e + 1 


63 


e + 1 


e 


e + 1 


e 


e + 1 


e 



(25) 



The first step has to be done following example 2. For all the rest the scheme of the first example can be applied. 
For the next 2e + 1 steps, we introduce vectors having three of the variables r, s, t, u positive and one strictly 
negative and r + s + t + u= —1, using the scheme of the first example: 



6 



a 


-1 


-2 


-k 


-(2e+l) 


bi 





1 


k-l 


2e 


b 2 


' 


o 


o ■•' 





b 3 















(26) 



We have constructed a set of Ad{d— 1) linearly independent vectors, which finishes the proof for the last remaining 
case. 

Thus, we have shown that the CGLMP-inequality satisfies not only condition 1, but also condition 2, and therefore 
it is a tight Bell inequality for all values of d. 



IV. COMPLETE CHARACTERISATION FOR GENERALIZED 3-OUTCOME CORRELATION 

FUNCTIONS 

In this section we give a way to simplify the data extracted from bipartite measurements, that is a generalization 
of the two-outcome correlation functions. With this simplification numerical work becomes possible and some results 
are obtained. 

Let us concentrate on the case d = 2, and let us label the two outcomes of all the observables by the numbers 1 
and —1, that is A a , B b = 1,-1. The correlation function of the variables A a and B b is defined as the average of their 
product, {A a B b ), for different realizations of the experiment. Now, instead of dealing with the 16-dimensional vectors 
with components 

P(A a = k 1 B b = s) k,s = 1,-1 a, b =1,2 (27) 

(like in section II and III), we could deal with the 4-dimensional vectors that have the correlation functions 

(A^), (A^), (A 2 B 1 ), (A 2 B 2 ) (28) 

as components. To transform a vector from the 16-dimensional space to its corresponding vector in the 4-dimensional 
space (28), we have to perform the projection 

(A a B b ) = P(A a = l,B b = l) + P(A a = -1, B b = -1) - P(A a = l,B b = -1) - P(A a = -1, B b = 1) , (29) 

in each of the four subspaces a,b = 1,2. It is known that the projection of a convex polytope is always a convex 
polytope [14]. Therefore, the set of vectors of correlation functions (28) achievable with local- realistic models is also 
a polytope, characterized by a new set of Bell inequalities. The CHSH-inequality [2] 

(AiBi) + (AiS 2 > + (A 2 B 1 ) - (A 2 B 2 ) < 2 , (30) 

is an instance. Even more, it was proven in [8], that all the non-trivial facets of this polytope are equivalent to 
the CHSH-inequality. Hence, the satisfiability of this inequality — and its equivalent ones — by a set of correlation 
functions is a necessary and sufficient condition for the existence of a local-realistic model giving these correlation 
functions. For the setting consisting of an arbitrary number of parties each possessing two dichotomic observables, 
all the Bell inequalities for correlation functions have also been obtained [10]. 

Dealing with correlation functions is much more simple than with joint probabilities, although in general, after 
doing the projection (29), some non-local-realistic joint probabilities can be projected into local-realistic correlation 
functions. It was proven in [8] that in the setting of the CHSH-inequality (d — 2) this do not happen, or in other words, 
correlation functions contain the same information that joint probabilities concerning the existence of a local-realistic 
model accounting for them. In what follows we show a generalization of the correlation functions that can be used 
with observables of more than two outcomes: Instead of using the 4eP-dimensional vectors of joint probabilities (1), 
let us deal with the 4<i-dimensional vectors with components: 

d-l 

P(A a -B b = n)=J2 p ( A a =n + j mod d, B b = j) . (31) 

j=o 

We say that the P(A a — B b = n)'s generalize the idea of correlation functions because they are also a projection 
of the joint probabilities and, for d = 2 (31) contains the same information than (29). Notice that this description 
is equivalent to say that only the information of the variables r,s,t,u (section III) is used. Because the CGLMP- 
inequality (11) can be written in terms of only r,s,t,u, the projection (31) leaves it unchanged, and this implies 
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that it is also a facet of the new polytope that yields the projection (31). This polytope is the convex hull of the d 3 
generators 

|4i - Bi> e |^i - b 2 ) e \a 2 - Bi> e \a 2 - b 2 ) a 1 ,a 2 ,b 1 ,b 2 = o,i...<i-i. (32) 

We have done an algorithm that finds all the facets of this polytope for any d. Contrary to what happens with the 
polytope of general joint probabilities, in this case the program runs quite fast for d = 3, and it turns out that, all 
non-trivial inequalities are equivalent to the CGLMP-inequality: 

Result: The satisfy ability of the CGLMP-inequality for d = 3 — and its equivalent ones — by a set of generalized 
correlation functions (31) is a necessary and sufficient condition for the existence of a local-realistic model accounting 
for them. 

Hence, we have a simple complete characterization of all local-realistic generalized correlation functions of three 
outcomes. For d = 4 the program takes too much time for an exhaustive search, but we have found that, there are at 
least three non-trivial inequalities which are not equivalent. This makes natural the assumption that for any d larger 
than three there is not a unique class of inequalities. 

V. CONCLUSIONS 

In this paper we have proven that the CGLMP-inequality is tight, by taking a geometric point of view. This is the 
first time that the tightness of a Bell inequality for an arbitrary number of outcomes is proven. There are not much 
criteria for evaluating the power of a Bell inequality, and tightness is a very objective one. 

We give a possible generalization of correlation functions that can deal with more than two-outcome experiments. 
With this definition a complete characterization of the polytope for d — 3 is numerically possible, and we have 
found that, all non-trivial facets are equivalent to the CGLMP-inequality. That is, for d = 3, the satisfyability of 
this inequality — and its equivalent ones — by a set of generalized correlation functions is a necessary and sufficient 
condition for the existence of a local-realistic model accounting for them. 

The fact that this inequality is not maximally violated by the maximally entangled state, for d > 3, could be due 
to the usage of only two Von Newman observables per party. Perhaps, increasing the number of observables per party 
or allowing for POVMs could remove this bizarre situation, although this is an open problem. 

The author is most grateful to J. I. Latorrc for its help in the design of the algorithm and many interesting 
discussions. The author also thanks A. Acfn, E. Jane and E. Rico for their useful comments and suggestions on 
the final manuscipt. This work is financially supported by the following projects: AEN99-0766, 1999SGR-00097, 
IST-1999-11053. G.V. and the grant 2002FI-00373 UB. 

VI. APPENDIX 

Here there are the two lemmas that have not been written in the main body of the paper: 

Lemma 1: Here, we prove that the non-signaling (2) and normalization (3) conditions give us a set of 4d linearly 
independent equations. 

To achieve this, we construct a set of Ad equations by adding one at every step which contains a coordinate that 
the previously introduced equations did not contain. First, we introduce the following equations: 

d-l d-l 

^2P(A 1 =k,B 1 = s) = ^2P(A 1 =k,B 2 = s) k = 0, 1, . . . d - 2, (33) 

taken from the set (3) — notice that the equation corresponding to k = d — 1 is not taken — . It is easy to see 
that for a particualr value k = ko the coordinate P{A\ = ko, B\ = d — 1) appears only in one of the equations 
— the one concerning A\ = fc — and therefore (33) is a set of d — 1 linearly independent equations. Second, we 
introduce the set of equations obtained by substituting A\ by A 2 in (33). As before, for every particular value 
k = k the coordinate P{A 2 = ko, B\ — d — 1) appears only in one of the equations introduced until now. Next, 
we do the same but with B\ and B 2 . After these operations, we have constructed a set of A(d — 1) linearly 
independent equations. Finally, we add the four normalization constrains (2). Notice that in each of these four 
equation a coordinate like P(A a = d — 1, B^ = d — 1) appears for the first time. This happens because we have 
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excluded the four equations corresponding to the value k = d — 1 from the non-signaling condition. It can be 
seen that these four equations are linear combinations of the ones in the set we have constructed. Hence we 
have proved that among the equations (2-3) there are only 4d which are linearly independent. 

Lemma 2: Here, we prove that the affine hull of the polytope has dimension 4d(d — 1). 

By construction all the generators satisfy the normalization (2) and the non-signaling (3) conditions, hence the 
whole polytope belong to the affine space found in Lemma 1. A generic generator with components SA a ,k °~B b ,s 
can be written in this way: 

G (Ai)A2>Bu b 2) = (L4i) 8 \A 2 )) ® (l-Bi) \B 2 )) . (34) 

Where the symbol \A) represents a <i-dimcnsional column vector with a 1 in the A th component and a in the 
rest. Let us concentrate for a while on the set of vectors 

(|Ai> © \A 2 )) E TZ 2d A 1 ,A 2 = 0,l...d-l. (35) 

It is easy to see that the following 2d — 1 vectors are linearly independent: 

(|o> e |o» , (|o)e|i», (|o)e|2» ... (|o)©|d-i», (|i>e|d-i», (|2)e|d-i» ... (\d- 1) © \d- 1». 

(36) 

Thus, we can say that in the set (35) there are at least 2d — 1 linearly independent vectors, or in other words, 
the linear span of this set has dimension larger or equal than 2d — I. The generators (34) are tensor products of 
vectors of the form (35). This implies that the linear span of the generators has dimension larger or equal than 
(2d — l) 2 . Because the null vector cannot be obtained by any affine combination [13] of the generators, their 
affine hull has one dimension less that their linear span. This means that the dimension, h, of this affine space 
is larger or equal than 4d(d — 1). But we know from Lemma 1 that the polytope belongs to an affine space of 
dimension 4d(d— 1), which implies that the dimension of its affine hull, h, cannot be larger than 4d(d— 1). Then, 
the only possibility is that the dimension of the affine space generated by the polytope is exactly h = 4d(d — 1), 
which concludes the proof. 
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